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Abstract 

We have studied the spin dynamics in Pro.63Sro.37Mn03 above and below 
the Curie temperature Tc = 301 K. Three distinct new features have been 
observed: a softening of the magnon dispersion at the zone boundary for 
T < Tq, significant broadening of the zone boundary magnons as T — > Tc, and 
no evidence for residual spin- wave like excitations just above Tc- The results 
are inconsistent with double exchange models that have been successfully 
applied to higher Tc samples, indicating an evolution of the spin system with 
decreasing Tc- 
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The revival in the study of manganites has led to a reexamination of the unique cou- 
pling between magnetism and charge transport in these materials. We focus on perovskite 
manganites with a transition from a high temperature paramagnetic insulator to a low tem- 
perature ferromagnetic metal at the Curie temperature Tc- The samples that exhibit this 
behavior have been partially hole doped away from a parent antiferromagnetic insulator, 
such as LaMnOs, by divalent substitution on the La site, such as Lao.7Cao.3Mn03 [jlj. The 
Mn 3d levels, split by the oxygen octahedral crystal field to a lower energy t 2 g triplet and a 
higher energy e g doublet, are filled according to Hund's rule such that all spins are aligned 
on a given site by a large intra-atomic exchange Jh- Electronic conduction arises from the 
hopping of an electron from Mn 3+ to Mn 4+ with electron transfer energy t. This results in 
the ferromagnetic double exchange interaction between localized S=3/2 spins (the core t 2g 
triplet) mediated by the hopping e g electron ]2]||. 

Recently, it has been shown that although the highest Tc materials (Tc > 325 K) are rea- 
sonably described by the double exchange model, with decreasing Tc (and reduced electronic 
bandwidth), the dramatic magnetotransport properties and increasingly first-order transi- 
tion require the incorporation of a strong Jahn- Teller based phonon coupling This 
has been corroborated by a number of studies showing the growing importance of static 
and dynamic lattice distortions as the zero temperature insulating state is approached [^,|| . 
Thus far, most studies of the spin dynamics have focussed on high Tc samples which have 
been shown to be consistent with the double exchange model PHTT]. With decreasing Tq 
however, important deviations may occur, as indicated by the development of a prominent 
diffusive central peak near T c in La .67Cao.33Mn0 3 (T c = 250 K) fll2"| . By studying the spin 
dynamics in a sample with reduced Tc, we can test if the spin system remains well described 
by simple double exchange (allowing for renormalizations of t and </#), or just as for the 
charge dynamics, the lattice coupling must be explicitly considered. 

In this paper we report a neutron scattering study of the magnetic dynamics in 
Pr 63 Sr 37 Mn0 3 above and below T c = 301 K, focussing on the zone boundary magnons. At 
10 K (0.03Tc), the magnons are well-defined for all wave- vectors q and exhibit a significant 
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deviation from the dispersion associated with the simplest local moment description. At 265 
K (0.9Tc), the dispersion relation has uniformly softened but maintains its q dependence, 
and significant broadening of the short wavelength magnons is observed. At 315 K (1.05Tc), 
just 14 K above Tc, there is no evidence for an inelastic magnetic peak at any q. We argue 
that these observations are at odds with double exchange calculations of the spin excita- 
tions (including leading order self-energy corrections), indicating qualitatively different spin 
dynamics in lower Tc samples. 

We used the floating zone method to grow a large single crystal of Pr .63Sr .37MnO3 with 
a mosaic spread of 0.5° full-width at half-maximum (FWHM). The neutron scattering mea- 
surements were carried out on the HB1 triple-axis spectrometer at the High-Flux Isotope 
Reactor, Oak Ridge National Laboratory. The collimations were, proceeding from the reac- 
tor to the detector, 50'-40'-S-40'-70', and the final neutron energy was fixed at Ef=13.5 or 
30.5 meV. The monochromator, analyzers and filters were all pyro lytic graphite. The sample 
was slightly orthorhombic and twinned at the measured temperatures. Even so, we assumed 
a cubic lattice (a Q = 3.86 A) because we could not resolve the effects of the twinning within 
the spectrometer resolution. 

Figure 1 shows the magnon dispersion along the [0,0,1], [1,1,0], and [1,1,1] directions at 
10 K, and along [0,0,1] at 265 K (0.9Tc). Figure 2 shows some of the constant q scans from 
which Figure 1 was obtained. The experimentally observed magnon peaks were accurately 
described by Gaussian fits. To verify that our results are not instrumental artifacts, we 
performed a 4-dimensional Monte Carlo convolution of the experimental resolution function 
1131 and the theoretical neutron scattering cross-section for damped spin waves. We verified 
that the peak positions were relatively unchanged, especially near the zone boundary saddle 
point, where significant resolution-induced shifts could occur. Our calculations place an 
upper bound of 0.5 meV for the shift in the peak position due to resolution effects. Another 
possible source of misinterpretation is that the broad peak at (0,0,1.4) at 265 K is actually 
a composite of phonon as well as magnon peaks. This was ruled out by measuring the same 
point in the next zone, which showed the same lineshape, down in intensity by the magnetic 
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form factor. 

Focussing first on the low temperature dispersion, a new feature we have observed is 
the significant softening at the zone boundary, seen in all directions. The Heisenberg spin 
Hamiltonian, H = — J2ij JijSi • Sj, couples the spins at site Ri and Rj by Jy. In the linear 
approximation, the spin wave dispersion relation is given by hu(q) = A + 2S(J(0) — J(q)), 
where J(q) = J2j exp[iq • (Ri — Rj)}- A allows for small anisotropies. The solid line in 
Figure 1 is the outcome of a fit for only nearest neighbor interactions for £ < 0.2 (the cubic 
zone boundary occurs at £ = 0.5), resulting in A = 1.3 ± 0.3 meV and 2SJ\ = 8.2 ± 0.5 
meV. Although this describes the data for £ < 0.2, the zone boundary magnons are missed 
by 15-30 meV. By comparison, the magnon dispersion in La .7Pb .3MnO 3 (To = 355 K) was 
found to be well described by only nearest neighbor interactions of similar strength |10| . 



The dashed lines in Figure 1 are a fit to the full data set including up to fourth neighbor 
interactions, resulting in A = 0.2 ± 0.3 meV, 2S J x = 5.58 ± 0.07 meV, 2SJ 2 = -0.36 ± 0.04 
meV, 2S J 3 = 0.36 ± 0.04 meV, and 2SJ A = 1.48 ± 0.10 meV. This accurately follows the 
dispersion except near the [0,0,1] zone boundary, which requires the next Fourier term in this 
direction, corresponding to J§. Nevertheless, the fourth nearest neighbor fit quantifies the 
remarkable result that additional extremely long range ferromagnetic couplings are required. 
Although J2 and J3 were necessary to fit the data, the more important correction to nearest 
neighbor coupling is J4. The long range and non-monotonic behavior of J(q) required by 
the data seems to rule out a simple Heisenberg Hamiltonian. 

On warming, the dispersion relation uniformly softens, as can be seen in Figures 1 
and 2 for the [0,0,1] branch at 265 K. At £ = 0.14, near the zone center, the magnon 
peak shows no substantial changes other than decreasing to lower energy as temperature 
is increased. In addition, it is resolution-limited at both 10 K and 265 K. Just above Tc, 
at 315 K, there is no evidence for the magnon peak as expected for the long- wavelength 
excitations. As the Brillouin zone is traversed (Figure 3), magnon lifetime effects become 
apparent. In particular, on approaching the zone boundary at £ = 0.5, the 10 K linewidth 
is substantially larger than the experimental resolution width. On warming to 265 K, the 



deconvolved widths, which we quote as full-widths at half-maximum throughout this paper 
and its figures, are nearly doubled, from their 10 K value of 8.4±0.5 meV to 13.2±1.9 meV. 
Again, at 315 K> Tc, no obvious magnon peak remains. 

An important issue to address is the role of the magnetic Pr ions, which display the 
crystal electric field (CEF) level shown in Figure 2. A recent powder neutron diffraction 
study of a manganite with similar Pr concentration observed a refined ferromagnetic moment 



of ~ 0.5/x.b at low temperature |14]1 . Thus the proper description of the total spin system 
involves a two component non-Bravais lattice. A bound on possible Pr-Mn interactions 
can be placed by examining the crossing of the CEF level and the magnon dispersion, for 
example near £ = 0.2 in the [0,0,1] branch. Even for weak coupling, this region is susceptible 
to the effects of mixing due to energy degeneracy. Experimentally, we find that the CEF 
level varies smoothly through the crossing of the magnon branch without dispersion. The 
integrated intensities of both the CEF excitation and the magnons show no anomaly at their 
crossing, and in particular, the observed magnon linewidth presented in Figure 3 shows no 
feature at the crossing. Also, the magnon dispersion passes smoothly through the CEF level 
just as in Lao.7Pbo.3Mn03, which does not contain a magnetic rare earth ion (Figure 1). 
Thus, our experiments show no evidence for Mn-Pr coupling with an effect on the Mn spin 
dynamics. 

To place our study of Pr .63Sr .37MnO3 in proper context, it is useful to compare with 
more conventional ferromagnets. Perhaps the most relevant comparison is to the europium 
chalcogenides, insulating ferromagnets that are considered ideal Heisenberg spin systems 
with extended exchange interactions (next nearest neighbor interactions are significant) 



15HT7|. In EuO, both J\ and J2 are ferromagnetic, whereas in EuS, J2 is a competing 
antiferromagnetic interaction with J\ (higher order terms are negligible). Although similar 
to Pr .63Sr .37MnO3 in this respect, the abrupt zone boundary softening we have observed 
(not seen in the europium chalcogenides) requires a large J4 while J2 and J3 are relatively 
small. Despite the large difference in T c for EuO (Tc = 69.2 K) and EuS (Tc = 16.6 K), 
the low frequency spin stiffness, magnon bandwidth, and Tc are self-consistently related via 



simple mean field theory. Comparison of the two manganite samples with Tc = 355 K and 
Tq = 301 K (Figure 1) shows that at low energy, the magnon dispersions are almost identical, 
while the zone boundary softening of the lower Tc sample results in a ~ 69% decrease in 
the magnon bandwidth from ~ 108 meV to ~ 74 meV. Although Lao.7Pbo.3Mn03 is roughly 
consistent with mean field theory, with decreasing Tc this breaks down because the spin 
stiffness does not diminish - the additional exchange constants have been exactly balanced 
by a reduction in the effective J\. 

The spin dynamics of both EuO and EuS have been extensively studied |L^,[n|, and 
the damping near Tc in both cases were observed to vary consistent with predominant 
magnon-magnon damping [ 20|,pT |, with the linewidth V varying as V oc g 4 ln 2 (/c£T/ftu; q ) for 
hujq <C k^T , and T oc g 3 for hujq ^> k#T . Although theoretically valid only in limiting 
extremes, and for (a a q) 2 <C 1, T(q) was observed to rise smoothly across hu^ = k B T in 
EuS. By contrast, we have observed the lack of a q dependence near the zone boundary for 
£ > 0.35. At present, it remains unclear whether our results indicate another important 
damping channel (quasiparticles, lattice distortions), or whether a revised spin wave theory 
which takes account of the unusual spin wave dispersion for Pr .63Sr .3 7 MnO3 is sufficient. 
Finally, in both EuO and EuS, at short wavelengths, a magnon-like peak was observed far 
above Tc, whereas we did not observe such a feature just above Tc in Pr .63Sr .37MnO3 for 
any £ in the [0,0,1] branch. 

We now discuss possible explanations of our results. Within a Stoner model, the conduc- 
tion band in the ferromagnetic state is exchange split into a majority and minority band. 
The magnon dispersion enters the Stoner continuum at finite q and u, where quasiparticle 
damping of spin waves occurs. The ferromagnetic ground state of the manganite, however, 
is rather unique in that there is complete separation of the majority and minority band 
by a large J#. Thus at low temperatures, the entire spin wave dispersion probably lies 
below the triplet electron-hole pair excitation continuum. At temperatures approaching Tc, 
however, the ordered moment is decreased and the carriers are no longer fully polarized, 
bringing triplet electron-hole pair excitations to energies within the spin wave band. These 



excitations can be the decay products of low energy magnons, and as their probability rises 
with temperature, give rise to temperature-dependent magnon lifetimes. 

The damping mechanism described above can be encapsulated in the imaginary part 
of the magnon self-energy. The real part contains the magnon dispersion, including its 
deviations from a simple cosine form. The magnon dispersion incorporating the lowest 
order self-energy correction has been calculated in a Kondo lattice model Hamiltonian with 



ferromagnetic spin-electron exchange 0]. For J#/t = oo, the result was a simple cosine 
dispersion consistent with a nearest-neighbor ferromagnetic Heisenberg exchange. For finite 
values of Jn/t, there are deviations from the simple result, but the deviations do not fit our 
experimental observations. In particular, the reduction of the low frequency spin stiffness 
is more significant than the reduction of the total bandwidth, contrary to the experimental 
results emphasized earlier. Beyond these inconsistencies, however, the most important point 
is that with decreasing Tc, t is reduced. Jh is presumably unchanged, being an intra-atomic 
energy. Thus with decreasing T c , corrections of order t/J H should decrease, not increase. 
Therefore, lowest order perturbative corrections in tj Jh do not capture our experimental 
results. 

The Kondo-lattice calculations discussed above have two limitations: the rate of con- 
vergence is not well understood, and the effects of varying the hole concentration x are 
not explored. These issues have been studied in one dimension by exact diagonalization 
of double exchange coupled spin rings for large Jh p2|j23[| . Deviations from a simple co- 
sine dispersion are predicted, particularly for extremely large or small doping (x < 0.2 and 
x > 0.7). For intermediate carrier concentrations (as studied here), however, a cosine band 
is recovered. For the two samples compared in Figure 1, the nominal carrier concentration 
is not significantly different, leading to the conclusion that a degree of freedom beyond x, t, 
or Jh is needed to account for our data. 

One much discussed parameter beyond the purely electronic parameters is the electron- 
lattice coupling, which is worth considering because of the large and strongly temperature- 
dependent mean square displacements (u 2 ) of Mn and O near Tc 0H- The bare overlap 



integral between the Mn spins is highly sensitive to the geometric arrangement of the Mn- 
O-Mn linkage, particularly due to the directional nature of the Mn d and p orbitals. The 
softening of the well-defined zone boundary magnons at low temperature would then be due 
to slow (on the scale of the zone boundary spin wave energies) fluctuations towards doubling 
of the nuclear unit cell. That fluctuations of this type might exist is clear from the myriad of 
lattice and charge-ordering instabilities in the phase diagrams of the manganites. Coupling to 
the lattice can also contribute to the magnon lifetime effects in the manganites. In particular, 
the rapid rise in (u 2 ) with temperature could give rise to a broadening distribution of effective 
exchange constants, leading to the pronounced broadening of the short wavelength magnons 
near Tq. Long wavelength magnons, of course, would not be affected as the variations in J 
would be averaged out. This scenario is qualitatively consistent with the T c dependence of 
the spin dynamics, in that the highest Tq samples do not show the strong feature in (u 2 ) 

We thank B. Batlogg, J. A. Fernandez-Baca, N. Furukawa, Y. B. Kim, P. B. Littlewood, 
J. W. Lynn, P. Majumdar, A. Millis, B. I. Shraiman, S. H. Simon, and C. M. Varma for 
useful discussions. HYH and GA would like to thank the HFIR staff at ORNL for their 
hospitality and assistance. The work at ORNL was supported by the U.S. DOE under 
contract No. DE-AC05-96OR22464 with Lockheed Martin Energy Research, Inc. 



S 



REFERENCES 

[1] G. H. Jonker and J. H. Van Santen, Physica 16, 337 (1950). 
[2] C. Zener, Phys. Rev. 82, 403 (1951). 

[3] P. W. Anderson and H. Hasegawa, Phys. Rev. 100, 675 (1955). 

[4] H. Y. Hwang, S-W. Cheong, P. G. Radaelli, M. Marezio, and B. Batlogg, Phys. Rev. 
Lett. 75, 914 (1995). 

[5] H. Roder, J. Zang, and A. R. Bishop, Phys. Rev. Lett. 76, 1356 (1996). 

[6] A. J. Millis, B. I .Shraiman, and R. Mueller, Phys. Rev. Lett. 77, 175 (1996). 

[7] P. Dai, J. Zhang, H. A. Mook, S.-H. Liou, P. A. Dowben, and E. W. Plummer, Phys. 
Rev. B 54, 3694 (1996). 

[8] P. G. Radaelli, M. Marezio, H. Y. Hwang, S-W. Cheong, and B. Batlogg, Phys. Rev. B 
54, 8992 (1996). 

[9] M. C. Martin, G. Shirane, Y. Endoh, K. Hirota, Y. Moritomo, and Y. Tokura, Phys. 
Rev. B 53, 14285 (1996). 

[10] T. G. Perring, G. Aeppli, S. M. Hayden, S. A. Carter, J. P. Remeika, and S-W. Cheong, 
Phys. Rev. Lett. 77, 711 (1996). 

[11] N. Furukawa, J. Phys. Soc. Jpn. 65, 1174 (1996). 

[12] J. W. Lynn, R. W. Erwin, J. A. Borchers, Q. Huang, A. Santoro, J-L. Peng, and Z. Y. 
Li, Phys. Rev. Lett. 76, 4046 (1996). 

[13] M. J. Cooper and R. Nathans, Acta Cryst. 23, 357 (1967). 

[14] P. G. Radaelli, G. Iannone, D. E. Cox, M. Marezio, H. Y. Hwang, and S-W. Cheong, 
Physica B, in press. 

[15] L. Passell, O. W. Dietrich, and J. Als-Nielsen, Phys. Rev. B 14, 4897 (1976). 

9 



[16] H. G. Bohn, W. Zinn, B. Dorner, and A. Kollmar, Phys. Rev. B 22, 5447 (1980). 
[17] H. A. Mook, Phys. Rev. Lett. 46, 508 (1981). 

[18] O. W. Dietrich, J. Als-Nielsen, and L. Passell, Phys. Rev. B 14, 4923 (1976). 
[19] H. G. Bohn, A. Kollmar, and W. Zinn, Phys. Rev. B 30, 6504 (1984). 
[20] V. G. Vaks, A. I. Larkin, and S. A. Pikin, Sov. Phys. JETP 26, 647 (1968). 
[21] A. B. Harris, Phys. Rev. 175, 674 (1968). 

[22] J. Zang, H. Roder, A. R. Bishop, and S. A. Trugman, J. Phys.: Cond. Mat. 9, L157 
(1997). 

[23] T. A. Kaplan and S. D. Mahanti, J. Phys.: Cond. Mat., in press. 



10 



FIGURES 

FIG. 1. Magnon dispersions for (0,0, £), (£,£,0), and (£,£,£) (where £ = 0.5 is the cubic zone 
boundary) at T = 10 K and 265 K. The solid line is a fit to a nearest-neighbor Hamiltonian for 
T = 10 K and £ < 0.2. The dashed line is a fit for all data including up to fourth nearest neighbors 
at T = 10 K. The dotted line is the corresponding four neighbor fit for T = 265 K. Also shown in 
squares are data for Lao.7Pbo.3Mn03 at 10 K (from Ref. 10). 

FIG. 2. Constant-q 1 scans at two different wave vectors with the left panel close to the zone 
center and the right panel close to the zone boundary. There is a dispersionless crystal electric 
field (CEF) level at ~12 meV from Pr (shown as dash-dot line). The solid circles are a constant-g 
scan of the CEF level at (0,0,1), the zone center. At 265 K, the intensity of the CEF level drops 
to undetectable levels and has therefore been ignored in the data analysis at this temperature. 
The T = 315 K data has been fitted with a simple Lorentzian line-shape (including a Gaussian). 
The other data have been fitted to a damped harmonic oscillator, including 4D convolution of the 
instrumental resolution as described in the text. The dashed line is the instrumental response to 
spin- waves with infinite lifetimes and the dispersion shown in Figure 1. 

FIG. 3. Magnon line widths in the [0,0,1] direction extracted from Gaussian fits with a sloping 
background. The solid and dashed lines are guides to the eye. 
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